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The g-state Potts model with long-range interactions that decay as 1/r" subjected to an uniform 
magnetic field on d-dimensional lattices is analized for different values of q in the nonextensive 
regime < a < d. We also consider the two dimensional antiferromagnetic Ising model with the 
same type of interactions. The mean field solution and Monte Carlo calculations for the equations of 
state for these models are compared. We show that, using a derived scaling which properly describes 
the nonextensive thermodynamic behaviour, both types of calculations show an excellent agreement 
in all the cases here considered, except for a = d. These results allow us to extend to nonextensive 
magnetic models with < a < d a previous conjecture which states that the mean field theory is 
P^ ' exact for the Ising one. 

C/3 






a 



o 



0^ 



X 



I. INTRODUCTION 



i-rt ! Microscopic pair interactions that decay slowly with the distance r between particles appear in different physical 

P3 ' systems. Typical examples are gravitational and Coulomb interactions, where the potential decays as 1/r. Several 

O " other important examples can be found in condensed matter, such as dipolar (both electric and magnetic) and 

1^1 - Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions, both proportional to l/r^. Effective interactions with a power 

I law decay l/r" , with some exponent a > 0, appear also in other related problems such as critical phenomena in highly 

^"^ ■ ionic systems |l), Casimir forces between inert uncharged particles immersed in a fluid near the critical point and 

J^ phase segregation in model alloys [^ . 

It is known that some of these systems can exhibit nonextensive thermodynamic behaviour (see Refs. ^Q and 
^-T) ] references therein). In other words, for small enough values of the ratio a/d the free energy F = —\nZ/(3, with 
\^ • Z = Trexp{—f3H) {H being the Hamiltonian of the system, d the dimensionality and f3 = l/fcsT), grows faster than 
^P I the number TV of microscopic elements when N — > cx), and the thermodynamic limit is not well defined. 

In a recent communication ||5|] two of us analized the thermodynamics associated with the long range (LR) ferro- 
magnetic Ising Hamiltonian 
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J(.nj) = — (J > 0; « > 0) (2) 



where Ty is the distance (in crystal units) between sites i and j, and where the sum ^,^ > runs over all distinct pairs 
of sites on a d-dimensional hypercubic lattice. It was shown |q| that the quantities per particle: free energy /, internal 
energy u, entropy s and magnetization ttt, of a finite model of N spins behave according to Tsallis conjecture [p[ for 
N ^ 1. These quantities present, in the presence of an external magnetic field h, the following asymptotic scaling 
behaviours: 

u{N,T,h) ^ N*u'{T/N*,h/N*) (3) 

f{N,T,h) ^ N* f{T/N*,h/N*) (4) 

s{N,T,h)r^ s'{T/N*,h/N*) (5) 

m{N,T,h) r^m'{T/N*,h/N*) (6) 



for all a > 0, where the functions u', /', s' and m' are the corresponding quantities associated with the same model but 
with rescaled coupling J'(rij) = J{rij)/N* (these functions are independent of the system size N) and the function 
N* (a) is defined as 

N*{a)^-^iN'-'^^'-l) (7) 

1 — a/a 

which behaves, for N ^ oo, as 

{— 7^^ for a/d > 1 

In TV for a/d =1 (8) 

__L^]^i-a/d foj. < a/d < 1 

For a > d the thermodynamic functions per site do not depend on N and the system is extensive (i.e., the thermo- 
dynamic limit exists). When < a < d, iV*(Q;) diverges for N ^ oo and the system is non-extensive. It was also 
presented numerical evidence |5| that for d = 1 the mean field (MF) theory becomes exact when < a < d. This led 
two of us to conjecture that the mean field theory might be exact for the nonextensive Ising model. 

In this work we extend the previous analysis and present new evidences of the exactness of the MF theory for 
< a < d for the one dimensional ferromagnetic g-state LR Potts model for different values of q (including the first 
order phase transition predicted by MF theory for q > 2) and also for the two- dimensional antiferromagnetic LR Ising 
model in an external field. We, thus, conjecture that the MF theory might be exact for any nonextensive magnetic 
model excluding the borderline case a = d, where there are probably corrections to the MF results (some previous 
experimental results W on dipolar ferromagnetic materials do not exclude this possibility). 

The outline of this paper is the following. In section II we analize the ferromagnetic q-state LR Potts model 
subjected to an uniform magnetic field h. First, we show that the previous analysis for the Ising model {q = 2) in the 
non-extensive region is straightforwardly extended to the generic q > 2 case. Then, in subsection II-A we derive the 
mean field solution of this model for arbitrary values of q, a and /i, in particular the MF predictions for the LR bond 
percolation which corresponds to the g — > 1 and /i — > 0"*" limit. In subsection II-B we compare the MF solution with 
our Monte Carlo simulation of the one-dimensional model for h — 0, q ~ 2,3 and 5 and different values of a. In section 
III we calculate the mean field solution of the two-dimensional Ising model with competing LR antiferromagnetic and 
short-range ferromagnetic interactions in an external field and compare them with the Monte Carlo results of Sampaio 
et al 1^]. Some comments and conclusions are presented in section IV. 

II. POTTS MODEL WITH LONG-RANGE FERROMAGNETIC INTERACTIONS 

In this section we address the LR ferromagnetic q state Potts model, i.e., we consider the Hamiltonian: 

H = -i ^J(ry)<5(a„a,)-/i^5(a„l) (a, = 1, 2, . . . , 9, V i) (9) 

i,j i 

where to each site, i, we associate a spin variable cr^, which can assume q integer values; the sum ^^ ■ runs over all 
distinct pairs of sites of a d-dimensional lattice of N sites (i j^ j); 6 is the Kronecker delta function, J{rij) is given by 
Eq.(0) and h is an external magnetic field in the a = 1 direction. The a — > cxd limit corresponds to the first-neighbor 
model. For g = 2 the a = limit corresponds, after a rescaling J -^ J/N, to the Curie- Weiss model. 

This model, in its plain formulation (a ^ cxd of Eq.(pf)) or in a more general one with many-body interactions, is at 
the heart of a complex network of relations between geometrical and/or thermal statistical models, like for example 
various types of percolation, vertex models, generalized resistor and diode network problems, classical spin models, 
etc (see [g[ and references therein). 

On the other hand, the Potts model with LR interactions has been much less studied. In the extensive regime 
a > d it presents a very rich thermodynamic behaviour, even in the one-dimensional case [|lO|-[I2|. To the best of our 
knowledge, no study has been carried out for the nonextensive regime < a < d. 



Let us introduce the sums (j)i{a) = J2i=jii ^/"''fj- ^ sufficient condition (and believed to be necessary [g3 14|) for the 
existence of the thermodynamic limit of this system is that 



Let us now consider a d-dimensional hypercube of side L + I and N = (L + l^, and let i = be the central site of 
the hypercube. We have that 

(f>{a) — lim <po{a). (11) 

When L ^ I {N ^ 1) 00 (q^) shows the following asymptotic behaviour B: 

0o(a)-^Q(«)2"7V*(a) (12) 

where N*{a) is given by Eq.(M) and Cd{a) is a continuous function B of a independent of N, with Cd(0) = 1 Vd. It 
can be proved that [|j: 

. . / 1 for < a < 1 , s 

^i(") = | -^C(«) fora>l • ^^^) 

where ^(z) is the Riemann Zeta function. From Eqs.([7[)-(|l2|) we see that the thermodynamic limit is well defined for 
a > d (where the system presents extensive behaviour), while for a < d the system becomes non-extensive. Following 
the same procedure as in Ref. [g| it can be shown that the scaling behaviours (|3|)-(y) of the thermodynamic functions 
hold Vg > 2 and Va > 0. For q = 2 the system undergoes a second order phase transition at finite temperature for all 
a > d when d> 2 |l^ and for 1 < a < 2 when rf = 1 |Q . For a ^ (i+ , the critical temperature shows the following 
asymptotic behaviour |15[| : 

kBT,^J(bia) (14) 

For d > 2 and a ^ d (short-range case ^) there exists a critical value qc such that the phase transition is a second 
order one when q < qdd) {qc = 4 for d = 2) and a first order one for q > qdd). For d = 1 and 1 < a < 2 Monte Carlo 
simulations U^ show that, for q > 2, there is a q-dependent threshold value adq) such that the transition is of first 
order when a < adq) and of second order above it. 

A. Mean field theory 

In order to develop a mean field version of Hamiltonian (0) we use Mittag and Stephen p^ ] spin representation for 
the Potts model, i.e., we associate to each site j a spin variable Xj which can take the values Aj = 1, u, uP' , . . . , w'~^, 
where w = e^'^^l'i is a gth root of unity. In other words, if the site j is in the state o then \,j = lo"^^ . Then, using 
the property 

g-i ^ A'^A'9-'==,5(A,A') (15) 

we can rewrite the Hamiltonian (O) as: 

^ = -^E^(-^.)E^'T'-^EE^' -^(^M (16) 

^ JJ l=\ ^ i 1=1 



where the constant term C{J^ h) is 



CiJ,h)^^J2j(^^^) + — ■ (17) 



The fraction of sites in the state a, Ug- = (1/-^) X^i "^ ^i'^i ~ ^) ^^ ^^ ^^^^ representation is given by 



1 

q 



l + |]c.«-'('^-i)(A') 
1=1 



(18) 



and the order parameter for a symmetry breaking in the a — \ direction is 



"^^T^' (1^) 



which can be written, using Eq.(nq), as: 

9-1 



m = — 

q 



^EW- (20) 

1=1 



The mean field solution for this model can be easily found from the Variational Method in Statistical Mechanics 
|18| by using a non-interacting trial Hamiltonian Hq given by 

N g-1 

^o = -''EE^' (21) 

i=l 1 = 1 

where 77 is the variational parameter to be found as a function of temperature. The variational free energy F is given 
by 

F = Fo + {H-Ho)o (22) 

= Fo + (^-:^)^^(A^)o-l5]j(r,,)|](A^Af' )o-C(J,/i) (23) 

^ i 1=1 ^ i,j 1=1 

where the free energy Fq associated with Hq is 

N 
Fo = --In [exp {pi^q - 1)) + (-7 - 1) exp {-pii)] (24) 

and (. . .)o denotes the canonical average using the Boltzmann measure proportional to exp (— /3_ffo)- 
Using equality (|l5|) one gets that 

(A,)o - (A2)o = . . . = (Ar')o - mo (Vi) (25) 

where the variational order parameter mg (defined by an equation similar to Eq.(EO)) is related to ry through: 



exp {(i-qq) - 1 

™o = ,0 ..f 7T (26) 

exp {fjriq) + [q - I) 

and from the property A''^' = (A')* it follows that 

(A^Af )„ = ml . (27) 

The minimization condition leads to 

T] = -[J(j){a)mo + h] (28) 

q 

which, combined with Eq.(pq) gives the following mean field equation for the order parameter toq: 

^ exp[P{J(j){a)mo + h)] - 1 
"^° exp/3[J0(a)mo + /i] + (g-l) ■ ^ ^ 

In the a — > cxD limit (short range interactions) we have 0(a) — > z, z being the coordination number of the lattice, 
and we recover, for h ~ 0, Mittag and Stephen's lO result. For q = 2 the Hamiltonian (||) is equivalent to the Ising 
one with long range interactions, providing that J^""" — 2J^^"'^ and /i-P°**« — 2h^^"^^. In this case Eq.(p9|) reduces 
to 

mo = tanh [/3/2( J (?!)(a)mo + h)] (30) 

and we recover the result from Ref. 0], i.e., the long range version of the Curie- Weiss equation which describes a 
second order phase transition for /i = at ksTc/ J''^'^"'^4'{a) = 1. 



The variational free energy F calculated at the minimum (condition (|28|)) gives the following mean field free energy 
per site /: 



/3/ = lng+^^^ + -|/3/i+-(<7-l)/3J0(a)m^ 

+ {q - l)l3hmn -[l + {q- l)mo] In [I + {q - l)mo] - {q - 1)(1 - mo) In (1 - mo)} (31) 



where we have used Eqs. (|23|)-(|29|). 

For g > 3 and h = the transition is of first order and it is easy to verify from Eq.(pl|) that at the critical 
temperature the order parameter m jumps from zero to the value (q — 2)/{q — 1). The critical temperature for q > 3 
is given by 

kTjJ=(l^) /(") , . (32) 

'' \q-lj 2\niq-l) > 

which recovers, in the a -^ cxd limit, Mittag and Stephen's result [Q7|. 

In the q ^r 1 limit, which corresponds to a bond percolation where the bond probability occupancy between any 
two sites i and j is vl^fji the order parameter probability percolation Poo (defined as the probability that a randomly 
chosen bond of an infinite lattice belongs to a cluster of infinite size) can be derived from [M 



Pooip) = 1 + Jim ^ (hm I- hm ^ln(Z)] j 
h^o+ dh L?— 1 oq AT^oo N J 



(33) 



where h — f3h, Z is the partition function of the Potts model with coupling constants J/rf, and p is the first neighbor 
bond probability given by p = 1 — exp {—J/kBT). 

One can easily show, from Eqs. (29), (pi|) and (pSh, that the probability percolation Poo(p) is, as expected [Oj, 



exactly the q ^ 1 and h ^ 0^ limit of the order parameter m-o, namely. 

Poo = mo(<7 ^ 1,7^ ^ 0+) = 1 - expl^^^P^moiq ^ l,h ^ 0+)] (34) 

Kb -I 

or, in terms of p, 

Pooip; «) = !-(!- p)-^(")^~(P'") (35) 

from which it follows that, due to the divergence of 0(a) in the non-extensive regime, Poojfi < p < 1;Q < a < d) = 1 
and, hence, the critical probability pc{0 < a < d) = in agreement with the exact result |14| . 

One can easily prove, from Eq.(|35|), that the percolation probability Poo of a finite system with N bonds presents 
an asymptotic scaling behavior similar to Eq.(|y), namely 

Poo{N,p) ^ P^ip*) {N»l) (36) 

where p* is the variable p calculated at T* = T/N*, namely, 

/ JN*\ 
P* = l-exp(-^) (37) 

and P^ is the probability percolation associated with the long range bond percolation whose bond probability oc- 
cupancies are p* /rij. Using this rescaled variable p* , the MF order parameter equation becomes, for A^ >> 1 (see 
Eq.(ll)): 

Pao{p*]a) = 1 - (1 _p*)C.(a)2°P^(p*;a) ^gg) 

which leads, for different values of a G [0, d], to monotonously increasing distinct order parameters as p* varies from 
the critical probabilites pl{a, d) = 1 — exp [— l/(Cd(a)2")] to p* = 1. However, if one introduces a more convenient 
variable, namely, 

.*^l-exp(-:^)=l-(l-/)^ (39) 



then all these MF probability percolation curves for different values of a and d coalesce into a single curve described 
by the equation 



Pooir*) = 1 ^ {1 ^ r*) 



*\P=o(r*) 



(40) 



The critical value r* where Poo{r*) vanishes is r* = 1 — exp (—1) = 0.63212 . . ., which leads to the MF critical 
probability 



Pcia) = 1 - exp 



-1 



0(a) 



Combining Eqs. (pi), (|T^) and (pi) one verifies that pda -^ d'^) vanishes asymptotically as: 

Pc{a -^ d+) 
which, in the particular case of d = 1, gives 



Cd(d)2'i^d ' 



(41) 



(42) 



p,{a^l+,d=l)^-{a-l) 
Notice that the asymptotic behavior (Eq. (Esf)) coincides with the lower bound for pc(l < ct < 2,d — 1) [Eol 



(43) 



B. Monte Carlo results 

We performed a Monte Carlo simulation using the heat bath algorithm on the one-dimensional Hamiltonian (Q) 
with h = and periodic boundary conditions for N = 300, 600 and 1200, for q = 2, 3 and 5, and different values of 
< a < 1. We calculated the magnetization per spin ( [lS| ) as a function of T* = T/N* for different system sizes and 
performed a numerical extrapolation for 1/iV — > 0. 

In Fig. 1 we compare the numerical results of the g = 2 case for m{T*) vs. 2kBT* /2°'J for different values of a 
with the mean field solution (equation (pO[)). We see that all the numerical curves fall into a single one in excelent 
agreement with the MF prediction. 
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FIG. 1. Monte Carlo extrapolated results (symbols) of 7n(T*) vs. 2fcsT*/(2" J) compared with the mean field solution (solid 
line) for q = 2 



In Fig. 2 we make the same comparison for the q — 3 and q — 5 cases. The sohd hnes represent the MF solution 
given by equation (p9) for h = 0. The dotted hnes in this figure indicate the mean field prediction for the critical 
temperature jumps ^3^). Again we observe, in both cases, an excelent agreement between our simulations and MF 
results for < a < 1, including the first order phase transition for q > 3. 
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FIG. 2. Monte Carlo extrapolated results (symbols) of m(T* 
line) for q = 3 and g = 5. 



vs. ksT* /{2°'J) compared with the mean field solution (solid 



III. ANTIFERROMAGNETIC ISING MODEL WITH LONG-RANGE INTERACTIONS 



We now consider an Ising model with competing LR antiferromagnetic and short-range ferromagnetic interactions 
in an external field, which is described by the Hamiltonian: 



h = -jf yi ^*^j + jY.^^^^^-^Y1 ^* (-^^ = ±1 V i), 



<ij> 



r 



(44) 



where J > 0, Ji? > and the smn X]<i j> runs over nearest-neighbor sites of a d-dimensional hypercubic lattice. The 
above Hamiltonian reduces, for J^ = 1 and J = 0.5, to the model studied by Sampaio et al [p| through Monte Carlo 
simulations. 

A mean field version of this model can be obtained by considering the Hamiltonian 



Hmf — — 2_^ ^Iff ^i 



(45) 



with 



'e// 



■J^Z^'^j+Jf H ™J + ^' 



r 



(46) 



J nn I 



where the sum ^ ■ ^^ ,- runs over all nearest-neighbor sites of i and 



Tr 



^MF 



{s.j{5,e-^^-(t^-»}. 



with 



Zai 



MF 



Tr 



{s.}^' 



-PHmf{{S.}) 



(47) 



(48) 



Now we consider the particular case of a square lattice. Dividing our lattice into two square interpenetrated 
sublattices A and B we can propose a solution of the form: 



m if i G A 
m^ if i € B 



Let us introduce the functions 



^^H")-Ei 



r 
jeA y 



(49) 



(50) 



'(«)-E;^ 



r 



ieA 



ieA 



with 



It can be easily seen that 



(f>^^Ha) + (f>^-^\a) = (j>ia) 



(/)(i)(a) = 2-"/2 0(a) 
(/.(2)(a) = (l-2-"/2)0(a). 

Then, substituting Eqs.(||) and (|5^) into Eq.(^) we obtain 

-J0(i)(a)m^ - J0(2)(q,)toB-|-4J^toB + /i if i e A 



(51) 



'^^ff \-J(j)'^^'>{a)m'^ ~J(t)^^\a)m^ +AJFm^ + h if ieB 



(52) 



Combining Eqs.(|52|), Eq.([45|) and Eq.(|47|) we arrive, after some algebra, to the following set of MF equations for the 
magnetization m — m + rn^ and for the staggered magnetization nig = m — m^ : 



sinh{2P [h- {J (j){a) -4 JF)m]) 



s/i(2/3 [/i-(J(^(a)-4JF)m])+cos/i(2/3 [(J0(a) (2i-"/2 -l)+4JF)m^]) 



(53) 



-sink (2/3 [(J (/.(a) (2i-"/2 - 1) + 4 JF)m,] ) 



^ cosh (213 [h-{J(l){a)-AJF)m]) + cosh{2P [(J(/.(a) (2i-"/2 - 1) + 4 7^)771^]) ' 
For /i 7^ it is easy to verify that the only solution of Eq. (|5J) and Eq. (^3|) is mg — and 

m = tanh [/? (h — (J (j){a) — 4 Jp)m)] 



(54) 



(55) 



In figure 3 we compare a numerical solution of Eq.(|5^) with the Monte Carlo data of Sampaio et al g] for a = 1, 
Jf = ^, J — 1/2, (3N*{a) — (0.3)~^ and different values of N. The function Cd{a) for d = 2 was evaluated 
numerically. We obtained, for a = 1, 6*2(1) = 0.8813 ± 0.0001. The comparison for other values of a < 2 gave similar 
results. 
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FIG. 3. Monte Carlo simulations of Sampaio et al. for the magnetization m vs. a rescaled magnetic field for square lattice 
sizes of 32 x 32 (circles), 48 x 48 (squares) and 64 x 64 (triangles) for a = 1, Jf = 1, J = 1/2 and T* = T/N* = 0.3. The solid 
line represents the MF solution given by Eq.(55) for N=64x64. 

A similar comparison for a = 2 is made in figure 4 (6*2(2) — 0.746 ± 0.001). The Monte Carlo data of Sampaio 
et al M do not agree very well with the MF magnetization, suggesting that corrections to the mean field result 
should be taken into consideration in this borderline case (where a = d). Notice that theoretical [pTl predictions and 
experimental results H obtained for the critical behavior of d = 3 uniaxial ferromagnets with exchange and strong 
dipolar interactions show that corrections to the MF behavior are needed in this a = d = 3 case. 

Summarizing this section, we verify that the Monte Carlo simulations for the equation of state of the d — 2 LR 
Ising antiferromagnet are in excellent agreement with the MF prediction when < a < 2. 
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FIG. 4. Monte Carlo simulations of Sampaio et al. for the magnetization m vs. a rescaled magnetic field for square lattice 
sizes of 32 x 32 (circles), 48 x 48 (squares) and 64 x 64 (triangles) for a = 2, Jf = 1, J = 1/2 and T* = T/N* = 0.3. The solid 
line represents the MF solution given by Eq.(55) for N=64x64. 

IV. CONCLUSIONS 

We have analyzed, in this paper, two long range spin models with power law decaying interactions {r~°') under an 
uniform magnetic field h: (i) the q-state LR Potts ferromagnet on d-dimensional hypercubic lattices (including the 
q ^ I and h ^ 0^ case of LR bond percolation) and (ii) the LR square Ising antiferromagnet with first-neighbor 
ferromagnetic interactions. Both models present non-extensive thermodynamic behaviors when < a < d, but their 
thermodynamic functions become finite when conveniently scaled variables are used Ig] . We have derived this scaling 
for g > 2 and have shown that the mean field probability percolation (i.e. the percolation order parameter) satisfies a 
similar scaling. The derived MF solution for the free energy of the LR Potts model led to spontaneous magnetization 
curves which agree very well with our Monte Carlo simulations for d — 1 and g = 2, 3 and 5 states for different values 
of < a < 1. An excellent agreement occurred also between the derived MF equation of state for the above d = 2 
antiferromagnetic LR model and the Monte Carlo simulations of Sampaio et al g for distinct values of < a < 2. 
Our results strongly suggest that the mean field theory is exact for LR non-extensive spin models v^rith 
< a < d. This conjecture is also supported by the fact that Bergersen et al [^ Monte Carlo simulations of the 
correlation function of the LR d=l Ising ferromagnet reproduce the MF result when < a < 1. Since our conjecture 
predicts, for the d = 1 LR Potts ferromagnet, a first order transition for q > 3 and < a < 1, it matches nicely with 
previous results |12] exhibiting first order transition for 1 < a < ac{q). 

Our results show that mean field behaviour is robust against variations of the range of interactions a within the 
non-extensive region, for a large class of magnetic systems. If our conjecture were true, this would have important 
practical implications: if you are considering systems with slow enough decaying interactions then you do not need 
sophsiticated approximations. 

It would be very interesting to extend the present analyis to more general non-extensive systems of interacting 
particles with long-range interactions. 
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